Phenomenological theory of superconductivity near domain walls in ferromagnets 
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We develop a phenomenological model of superconductivity near a domain wall in a ferromag- 
net. In addition to the electromagnetic interaction of the order parameter with the ferromagnetic 
magnetization, we take into account the possibility of a local enhancement or suppression of super- 
conducting pairing in the vicinity of the wall, and also a non-perfect transparency of the wall to 
electrons. It is found that the critical temperature of superconductivity near the domain wall might 
be substantially higher than in the bulk. 

PACS numbers: 74.25.Dw, 74.25. Op, 74.20.-z 



I. INTRODUCTION 

The recent discovery of superconductivity co-existing 
with ferromagnetism in UGe2 Q , ZrZn 2 Q , and URhGc 
have revived the long-standing interest to the inter- 
play of the two phenomena. Although the microscopic 
mechanism responsible for superconducting pairing in the 
presence of a ferromagnetic background is still subject to 
considerable controversy, some useful information about 
the properties of these systems can be obtained from a 
phenomenological approach, based on general symmetry 
arguments and the Ginzburg-Landau (GL) equations, see 
e.g. Ref. Q. A symmetry-based description of ferromag- 
netic (FM) superconductors has been developed in Refs. 

In contrast to the non-magnetic case, a phenomeno- 
logical theory of superconductivity in the presence of fer- 
romagnetism should include the effects of internal mag- 
netic induction B created by the macroscopic FM mag- 
netization M. That the coupling of the orbital motion 
of electron charges with B would suppress supercon- 
ductivity even in zero external field was pointed out by 
Ginzburg [T^ |. who argued that superconductivity could 
not exist in any of the FM materials known at that time. 
While valid in a single-domain ferromagnet, his argument 
should be taken with care if B is non-uniform. In order 
to minimize the magnetic field energy, it is energetically 
favorable for a macroscopic ferromagnet to break up into 
domains of size L with opposite directions of the magneti- 
zation M , separated by domain walls (DWs) of thickness 
I. Due to the spatial variation of the internal field, it is 
expected that the suppression of superconductivity in the 
vicinity of a DW is not as strong as in the bulk, as was 
first hypothesized in Ref. |l3| . The quantitative analysis 
shows 14] that the superconducting critical temperature 
near a DW is indeed higher than the bulk critical tem- 
perature deep inside the domain. The reason is that the 
effective GL potential near the wall has the same double- 
well shape as near a superconductor-insulator interface, 
leading to a local enhancement of superconductivity simi- 
lar to the surface sheath effe ct llq . The same model was 
recently extended in Ref. [l6j to calculate the critical 
temperature of the superconducting DW in an external 
magnetic field. 



In all previous treatments of the superconducting DW 
problem it was assumed that the superconducting cou- 
pling constant is the same throughout the system. How- 
ever, it can be expected that, whatever the microscopic 
pairing mechanism, the conditions for the formation of 
the Cooper pairs near the DW and far from it are differ- 
ent. The FM magnetization affects the electrons via the 
orbital interaction of the electron charges with the induc- 
tion B, the exchange band splitting, and also the Zeeman 
coupling of the electron spins with B. Because of the spa- 
tial variation of both the exchange field and the internal 
induction, all three mechanisms could change both the 
single-electron spectrum and the electron-electron inter- 
actions in the vicinity of the DW, compared to the bulk 
of the domains. In particular, the strength of the Cooper 
pairing could be locally enhanced or suppressed near the 
wall. Also, the DW creates a potential barrier for elec- 
trons travelling from one domain to the other. The possi- 
bility of a non-perfect transparency of the barrier should 
therefore be included in the theory. 

We assume that the thickness of the wall is small 
compared with the length scales of superconductivity: 
I <C £o "C Xl, where £o is the coherence length and Al is 
the London penetration depth. Then, in addition to the 
usual (i.e. single-domain) terms, the general phenomeno- 
logical GL functional includes interface contributions lo- 
calized near the DW, which describe the variation of the 
pairing strength and a finite transparency of the DW for 
electrons. These additional terms affect the boundary 
conditions for the superconducting order parameter at 
the wall |l7j . and the problem then becomes formally 
similar to that of the localized superconductivity near 
a planar crystalline defect, which has been extensively 
studied in the context of the twinning plane supcrcon- 
ductivityin high-T c cuprates and other superconductors 

The purpose of the present article is to study the phase 
diagram of the localized DW superconductivity in a gen- 
eral phenomenological model. The article is organized as 
follows: In Sec. [HI we formulate our model and derive the 
GL equations with appropriate boundary conditions. In 
Sec. IIIII we calculate the critical temperature of the DW 
superconductivity as a function of an external field. The 
evolution of the shape of the superconducting nucleus in 
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the non-linear regime below the critical temperature is 
studied in Sec. II VI 

II. THE MODEL 

We assume that the Curie temperature Tfm is higher 
that the critical temperature T$c and treat the system 
in the mean-field approximation, so the superconductiv- 
ity appears in the presence of a static FM background. 
According to Refs. 0, IE El El, the order parameter 
in a FM superconductor transforms according to one of 
the irreducible co-representations of the magnetic sym- 
metry group in the normal state. Here we assume that 
the order parameter has one component ip(r), which 
is always the case in, e.g., ZrZn2, assuming that only 
one of the exchange-split bands is superconducting (in 
principle, the pair scattering can induce order param- 
eters of the same symmetry on different sheets of the 
Fermi surface, leading to a more complicated form of the 
free energy functional than the one used below). The 
spin quantization axis is chosen parallel to M(r), so 
that the pairing in both domains occurs in, say, spin- 
up channel. Thus the superconducting order parameters 
are A + (k,r) = A^\k,r) = ip + (r)(j)(k), A_(fc,r) = 

(k,r) — T/j~(r)(f)(k), where 4>(k) is the basis func- 
tion of the chosen one-dimensional co-representation. 

The bulk critical temperature is the same on both sides 
of the wall and equal to T c q. We focus here on the case 
I <SC £o <C L, where £ is the superconducting coherence 
length, which allows us to consider a single wall separat- 
ing two domains with opposite directions of magnetiza- 
tion. The GL free energy density can be represented as 
a sum of the bulk contributions from the two domains, 
and the interface contribution due to the presence of the 
DW: F = F + + F-+F DW . Choosing the DW to coincide 
with the plane x = 0, we have 

F ± =a(T-T ca )\iij ± \ 2 + K\DiP ± \ 2 + ±[3\iP ± \ 4 , (1) 

where = ip(x > 0) and ijj— = t])(x < 0) are the order 
parameters in the right and left domains respectively, 
A = -iVj + (2-K/§ )Ai{r), V x A = B, $ = nhc/e is 
the magnetic flux quantum (the electron charge is equal 
to — e). The internal magnetic induction is given by B — 
4nM + if, where H is an external field. Since the DW 
thickness is much smaller that £o, we can approximate 
the non-uniform magnetization by a step-like function: 
M(r) — (0, 0, Mo signx). If the external field is directed 
along the z axis, we have A(r) = (0, -BqI^I + Hx,0), 
where ESq = \-kM§ (we assume Mo > 0). For simplicity, 
we neglect the anisotropy of the effective mass tensor in 
the gradient terms. 

The DW contribution, which is localized near x = 0, 
can be written in the following general form: 

F DW = [ 7 i(i^+i 2 + iV'-i 2 )-72«>;^+cv' + ) 

+<ts(V;V- - 1>-1>+)Mx), (2) 



with real coefficients 7$. We keep only the terms 
quadratic in tjj±, which are needed to obtain linear 
boundary conditions at x = 0, see below. The expres- 
sion J2J is consistent with the symmetry of the system: 
Fdw has to be real and invariant with respect to time 
reversal (ip± — > il>±) accompanied by the interchange of 
the two domains (+ <-> — ). The parameters of the GL 
functional depend on the domain magnetization M . In 
particular, in the limit M — » (no domain wall), all 

71 should vanish. While the first two terms in Eq. 

are similar to those discussed in the context of the local- 
ized superconductivity at planar crystalline defects [19(, 
the last one is unique to the present problem because it 
requires time-reversal symmetry breaking. It turns out 
however that we can choose 73 = without any loss of 
generality, because the second and third terms in Eq. 
can be combined together by rotating the order param- 
eter phases separately in each domain: ijj± — > e TlS / 2 ip±, 
where 8 = axg^s/^)- For the same reason, the sign of 

72 is not important, in particular one can always choose 
it to be positive: if 72 < 0, then the gauge transformation 
-0+ — > ip+, ip- — * e t7r i[)- changes 72 — > —72- 

In our phcnomcnological model we allow for the or- 
der parameter to be discontinuous across the domain 
wall, the rationale being that the DW creates a poten- 
tial barrier for electrons and therefore may behave sim- 
ilar to an extended Josephson junction. The interface 
terms |J5J then translate into the linear boundary con- 
ditions for tp±, see e.g. Ref. 0. The magnitude of 
the coupling between the domains, described by the 72 
term, depends on the transmission coefficient of the DW 
for electrons, which in turn is determined by the thick- 
ness and the structure of the wall. In the limit of zero 
transparency, 72 = and the domains are completely 
decoupled. The continuity of the order parameter is re- 
stored in the limit 71 = 72 — > +00, when the inter- 
face terms can be written as ~ 4>-\ 2 , which im- 
poses the condition ip+(x = 0) = tp—(x = 0) [and also 
ip' + {x = 0) = ip'_(x = 0), see below]. If 71 = 72 + ^7, with 
71,2 - ► +00 and a finite S^/, then the order parameter is 
continuous, and Fjjw = 6^\ip\ 2 6(x). Depending on the 
sign of 57, this term describes either a local enhancement 
(<$7 < 0) or suppression (Sj > 0) of superconductivity 
near the DW. 

Far from the wall, the magnetic induction is uniform 
and the field dependence of the superconducting critical 
temperature is given by the standard upper critical field 
expression T Ct t, u ik — T c q — 2ttK Bq/&ocl (one should keep 
in mind that T c q itself can be a function of the magneti- 
zation). Thus the superconductivity in the bulk is sup- 
pressed by the internal induction Bq even in the absence 
of the external field At T < T c ,b u ik, both domains 
are in the mixed state. 

Let us now calculate T c in the vicinity of the wall. 
We assume that the phase transition into the supercon- 
ducting state is of second order at all fields. To find the 
phase diagram in the external field H, we use the lin- 
earized GL equations obtained from Eq. |JI]|. with the 
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boundary conditions supplied by the interface terms J3J • 
The order parameter has the form ip±(r) — e lqv f±(x), 
where the functions f±(x) satisfy the equations 



K 



d 2 .f± 
dx 2 



K 



2n 



(H±B )x + q 



f± 



-a(T - T c0 )f ± = 0, 



(3) 



supplemented by the boundary conditions at x 

72/-, 



= 0: 



K ~T= 7i/-+ 
ax 



K^r- = 72/+ - 71/- 
ax 



(4) 



To simplify the notations, it is convenient to use dimen- 
sionless variables: x — x/Ib, Q = qls, and % — ^Ib/K, 
where Is — ■J $>o/2ttBq. Omitting the tildas, the GL 
equations take the form 

-f£ + [(l±h)x±q] 2 f ± + (r-l)f± = 0, (5) 

where h = H/Bq > is the dimensionless external mag- 
netic field, and r = l+a(T— T c o)$o/2ttBoK is the dimen- 
sionless temperature. It is easy to see that the zero-field 
bulk critical temperature T Ct b u ik corresponds to r = 0, 
while T c o corresponds to r = 1. The boundary condi- 
tions 10} become 



/iT7i/±±72/ T = 0. 



(6) 



As seen from Eqs. (|5j), the effective GL potential V{x) 
is described by two parabolas matched at x = 0. We have 
to find the absolute minimum of the ground state energy 
1 — r in this potential as a function of the parameter q. 
The solution of Eqs. (J5J that does not diverge (in fact, 
vanishes) at |x| — * 00, has the form 



f±(x) =C ± e~ x i/ 2 H„ ± (±X ± ), 



(7) 



where C 
tions 



m 



are some constants, H v (z) is the Hermite func- 
and 



V±= 2 



1 - r 

\v±h\ 



- 1 



X ± = y/\l±h\ X ± 



1 



l±h 



These expressions are singular at h = ±1, which re- 
quires special consideration, see below. Substituting 
f±(x) in the boundary conditions Jfjjl and using the fact 
that H' v {z) = 2vH v ^\(z), we obtain the solutions of two 
types. First, there always exists a trivial bulk solution 
located deep inside the domain with lower induction, for 
which q = +00 or —00 and 1 — r = |1± h\, thus the crit- 
ical temperature is simply T Ci buik{h) = 1 — min |1 ± h\ = 
1 — |1 — h\ (at h > 0). Second, we can also have a non- 
trivial solution localized near the domain wall, which is 
sensitive to the values of 7$ . The corresponding equation 
for r(h, q) is 



V (^) 



7i 



L v . q 



7i 



(8) 



Here we introduced L u {z) = 2vH v ^x{z) / H v {z) — z. The 
critical temperature T c (h) is given by the absolute max- 
imum of r(h,q) with respect to q at given h. Clearly 
T c (—h) = T c (h) [this is consistent with Eq. JHJ being 
invariant under h — > —h], therefore we assume h > 0. 

As mentioned above, the case h = 1 requires special 
care, because of the singularities in Eq. JSJl. Physically, 
this corresponds to a complete compensation of the inter- 
nal induction by the external field in the domain x < 0. 
While f + (x) in this case is still expressed in terms of the 
Hermite functions, f—(x) is given by a simple exponen- 
tial. The equation for r(h = l,q) is 



\Jt - 1 + q 2 



■ 71 



0, (9) 



where v = — (r + l)/4 [note that the critical temper- 
ature of the DW superconductivity is never below the 
bulk critical temperature, in particular r c (h = !)>!]. 



III. PHASE DIAGRAM 

A. Transparent DW 

Before discussing the solution of Eq. JSJ for arbitrary 
values of 71 and 72, let us first look at some important 
particular cases. If the domain wall is completely trans- 
parent for the electrons, then 71 = 72 + (57, 71,2 — ■> +00, 
and f + = /_. In this limit, Eq. © reduces to 



1 - h 



-y/\l-h\L u _ 



= 2 57. 



(10) 



In the simplest case Sj — 0, corresponding to the absence 
of the local enhancement or suppression of superconduc- 
tivity, the results of Ref. [16J are recovered from Eq. 
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pO|l. At zero external field we have v± = v = — r/2, and 
Eq. (|10|l becomes L v (q) — 0, whose numerical solution 
shows that the ground state energy 1 — r = 2v + l reaches 
its minimum at v = v* ~ —0.20 and q — ~ —0.77. In 
terms of the critical temperature, this gives r c (ft = 0) = 
-2z/* ~ 0.41 and T c ~ T c0 - 0.59(27rA' J B /$ a) > r c , wfc . 
Thus, superconductivity can exist in the vicinity of the 
DW at a temperature higher than T c ,t u ik- This effect has 
the same origin as the surface critical field H c3 At 
ft 5^ 0, the critical temperature first grows as a function 
of ft,, because the external field partially compensates the 
internal induction in one of the domains, and then finally 
starts to decrease at ft > 1, see Fig. 1. 
At ft = but 67 ^ 0, we have 

L-r/2(q) = h (11) 

If the superconductivity is locally suppressed, i.e. £7 > 
0, then the critical temperature is always lower than at 
6"/ = 0. In contrast, if the superconductivity is locally 
enhanced, i.e. #7 < 0, then t c is always higher than at 
<5 7 = 0. 

The results of numerical solution of Eq. ijTUjl are pre- 
sented in Fig. 1. The bulk critical temperature T c ,bulk (h) 
serves as the lower boundary for r c (ft) in the limit of 
strong suppression of superconductivity near the DW. In 
the opposite limit of large negative £7 (i.e. strong super- 
conductivity enhancement), the influence of the variation 
of the internal induction is less pronounced, so that the 
phase diagram becomes similar to that of a supercon- 
ducting plane defect in a non-magnetic crystal [Tg| - The 
order parameter is continuous and localized in the vicin- 
ity of the DW, see Figs. 2 and 3. 

Another notable feature of the phase diagram is that, 
in contrast to the upper critical field in the bulk, r c (ft) for 
the DW superconductivity is quadratic in field at ft — ► 0: 
r c (ft) — r c (ft = 0) = 0(h 2 ). This happens because the 
localization of the order parameter near the DW makes it 
possible to treat ft as a small perturbation, which cannot 
be done in the bulk case. 



B. Decoupled domains 

The opposite case 71 = 72 = corresponds to com- 
pletely decoupled domains with the order parameters 
obeying /±(0) = 0. As follows from Eq. ©, the crit- 
ical temperature in this case is determined by the equa- 
tion L v _ [(Zyjl — ft|/(l — ft)] = at ft > 0, which gives 
Tc (ft) = l-(2i/* + l)|l-ft| ~ l-0.59|l-ft|, so we have 
an iJ C 3-type behavior at all fields. 

C. General case 

Let us now find the critical temperature and the shape 
of the superconducting nucleus at zero external field for 
arbitrary 71 and 72. At ft = Eq. JBJ yields two equa- 
tions for the critical temperature: L_ T / 2 (q) = 7±, where 



7± = 71 ± 72. Recalling that 72 > 0, the highest criti- 
cal temperature is achieved when the right-hand side is 
7_ = £7 = 71 — 72. At Sj — 0, the transparent DW case 
is recovered, see Sec. IIII Al 

At arbitrary ft, 71, and 72, Eq. (jSJ) is solved numer- 
ically. A representative phase diagram is shown in Fig. 
4. In contrast to the transparent case, the order param- 
eter is now discontinuous across the DW, but it still is 
localized in its vicinity, see Figs. 5,6. The qualitative 
features of the phase diagram are quite similar to those 
discussed in Sec. IIII Al In particular, the crossover from 
the bulk regime to the strongly-localized superconduct- 
ing DW regime is controlled by #7 = 71 — 72 ■ 

IV. NON-LINEAR EFFECTS 

We have also studied the evolution of the shape of the 
superconducting nucleus as a function of temperature at 
zero external field, assuming a transparent DW with no 
local enhancement or suppression of superconductivity. 

At 71 = 72 = +00, the order parameter and its deriva- 
tive are continuous across the wall. Writing the order pa- 
rameter in the dimensionless form: ip(r) — ipoe lqv f (x) , 
where ipo = y/K/ (31%, the free energy becomes 

lr (I)* + (T - 1,/2+(|I|+ 9)2/2 + 5' 4 < 12 > 

(Fq = K 2 /f3lg), which yields the following non- linear GL 
equation: 

-f" + (\x\+q) 2 f + (T-l)f + f = (13) 

The solution of the linearized version of this equation 
in Sec. IIII Al gave the critical temperature r c ~ 0.41, 
corresponding to q« ~ —0.77. 

The non-linear equation l|13|) is solved numerically at 
< t < t c by the "shooting" method We start at 
some large value of |x| = xq 3> 1 on both sides of the DW 
and then integrate the equation towards a matching point 
near the origin (although in the absence of an external 
field, the solution is symmetric about the DW, the match- 
ing point was not necessarily chosen to be at x = 0). For 
each t and q, the starting value of /(±xo) = was cho- 
sen, while the starting slope of /'(±xo) was allowed to 
vary to find the solution, which is smooth everywhere. 
The upper and lower bounds on the initial slopes were 
estimated using the fact that at large \x\ the non-linear 
term in Eq. i|13[l is negligible, and f(x) decays approx- 
imately as a Gaussian. Then the bisection method was 
used to converge on the starting slope that guaranteed 
continuity of the wave function. At the next step, for 
each wave function an evaluation of the GL free energy 
(|12fl was made, and the value of q was chosen to minimize 
F. The final results were found to be independent of the 
choice of the matching point. 

Our results for the order parameter are presented in 
Fig. 6. The height of the wave function, i.e. f(x = 0), 
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grows approximately as \Jr c (h = 0) — r. As a check of 
the validity of our procedure, the temperature depen- 
dence of the parameter q was found to converge perfectly 
to the previously calculated value as r approaches 

T c (h = 0). 

V. CONCLUSIONS 

We have studied the phase diagram of a superconduct- 
ing domain wall in a ferromagnet using a phenomenolog- 
ical model featuring additional interface terms in the GL 
energy. These terms are introduced in order to account 
for the possible variation of the superconducting coupling 



constant near the DW, and also for a finite transparency 
of the wall to electrons. The critical temperature for the 
localized DW superconductivity is found to be always 
higher than in the bulk. Depending on the relation be- 
tween the phenomenological parameters in the interface 
terms, one can get an enhancement of the critical tem- 
perature above its value for a transparent wall, and also a 
non-monotonous dependence of the critical temperature 
on the external field. 
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FIG. 1: The superconducting phase transition lines for a 

transparent DW with a local enhancement or suppression of FIG - 2; The superconducting order parameter and the GL 
superconductivity, i.e. 71 = 72 + ^7, 71,2 -> +00. The dashed effective potential V(x) for h = 0.3, at 72 = 00, Sj = -0.5. 
line corresponds to the bulk critical temperature T Ct bulk(h) = 
1 — |1 — h\, the thick solid line - to a transparent DW case 
with 57 = 0, and the thin solid lines correspond to £7 = 0.5, 
0.3, -0.5, -0.7, -0.9, -1.1 (from left to right). 





FIG. 4: The superconducting phase transition lines for 71 = 
72 + #7, 72 = 0.2. The dashed line corresponds to the bulk 
critical temperature, the thick solid line - to a transparent 
DW case with £7 = 0, and the thin solid lines correspond to 
Sj = 0.5, 0, -0.1, -0.3, -0.5, -0.7, -0.9 (from left to right). 





FIG. 6: The superconducting order parameter and the GL 
effective potential V(x) for h = 1.2, at 71 = —0.3, 72 = 0.2. 




FIG. 7: The superconducting order parameter in the non- 
linear regime below r c , in the transparent DW case with no 
superconductivity enhancement [r c (h = 0) ~ 0.41)]. The 
curves correspond to r = 0.40, 0.30, 0.15, and 0.05 (from 
bottom to top). 



